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We obtain a model-independent expression for the Dalitz plot of semileptonic 
decays of polarized hyperons including radiative corrections to order a and 
neglecting terms of order aq/irMi, where q is the four-momentum transfer 
and M\ is the mass of the decaying hyperon. We specialize our results to 
exhibit the correlation between the charged-lepton momentum and the spin 
of the decaying hyperon. We present results for the three-body region of the 
Dalitz plot and for the complete Dalitz plot (which includes the four-body 
region). From these results we also obtain the corresponding radiative correc- 
tions to the integrated lepton spin-asymmetry coefficient. Our formulas are 
valid for charged as well as for neutral decaying hyperons and are appropri- 
ate for model-independent experimental analysis whether the real photon is 
discriminated or not. 
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I. INTRODUCTION 



The form factors of hyperon semileptonic decays, A — > Bivg, contain important informa- 
tion about the low-energy strong interactions of spin-1/2 baryons (A and B are such baryons 
and i and are the accompanying charged lepton and neutrino). Their experimental deter- 
mination requires the use of accurate formulas in the analysis of the measurements of several 
observables. An important one of these observables is the angular correlation between the 
spin §! of A and the direction 1 of the momentum of t. It is the purpose of this paper to 
calculate the radiative corrections (RC) to the Dalitz plot (DP) with the spin correlation 
§! • 1 exhibited explicitly. 

We shall obtain expressions that are suitable for a model-independent experimental anal- 
ysis. The model dependence of the virtual RC is handled following the approach of Sirlin 
IHII to the RC of neutron beta decay whereas the model dependence of the bremsstrahlung 
is controlled by the theorem of Low [|2| . In previous work we have discussed the unpolarized 
DP and the DP with the Si ■ p2 spin correlation kept explicitly Q (p2 is the direction of 
the momentum of the emitted baryon B) . It is not possible to derive the result for the spin 
correlation §i • 1 from the final result of Ref. [ffl, because all kinematical integrations, except 
for the i and B energies E and E2, respectively, were already performed. However, since we 
are going to follow the same approach of this reference, much of the work has already been 
advanced. 

The bremsstrahlung RC is a four-body decay whose DP covers entirely the DP of the 
three-body decay A — > Blvg. We shall refer to the latter as the three-body region (TBR) and 
to the non-overlap of the former and the latter as the four-body region (FBR). Even when 
no experimental arrangement has been made to detect and discriminate real photons, it is 
possible to eliminate the photons that belong to the FBR by energy-momentum conservation. 
Therefore, in calculating the bremsstrahlung RC we shall keep a clear distinction between 
these two regions. 

We shall also obtain the radiative corrections to the integrated lepton spin-asymmetry 
coefficient ai. As we shall see the distinction between the TBR and the FBR leads to a 
perceptible change in the RC to this asymmetry coefficient. 

Our results will be presented in two final forms. One where the triple integration over 
the real photon three-momentum k is left indicated and ready to be performed numerically. 
And another one, an analytical form, where such a triple integration has been performed. 
Both forms can be used to numerically cross-check on one another. However, the analytical 
result, although tedious to feed into a Monte Carlo program, leads to a considerable saving 
of computer time, because the triple integration does not have to be performed within the 
Monte Carlo every time E and E 2 or the form factors are changed. 

For the use of our results it is important that this paper be as self-contained as possible. 
In Sec. II we introduce our notation and conventions and we review the virtual RC; also the 
infrared divergence of this part is clearly separated. In Sec. Ill the real photon emission is 
calculated and separated into the contributions from the TBR and FBR. The calculation of 
Ref. H is adapted to the present case. The infrared divergence is extracted following Ref. || 
and its cancellation with the one of Sec. II is discussed in detail. Our first main result is 
established, allowing for the elimination (or not) of real photons from the experimental 
analysis. In Sec. IV we proceed to the analytical evaluation of the triple integration over 
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the bremsstrahlung three-momentum. Our second main result is established, also allowing 
for the experimental discrimination (or not) of real photons. In Sec. V we use the analytical 
result to obtain the RC to the asymmetry coefficient ag. In Sec. VI we make numerical 
evaluations for several HSD and also for the A+ — > Ae + u decay. We compare with other 
results available in the literature. Section VII is reserved to discuss and summarize our 
results. To make this paper self-contained we introduce three appendices. In Appendix 
A we give the amplitudes for virtual and bremsstrahlung RC, emphasizing how the model 
dependence is kept under control. In Appendix B we give the analytical expressions of all 
the coefficients, both new and of Refs. || and || required to compute the RC to the Sj ■ 1 
correlation. Finally, in Appendix C we review briefly the procedure of Ref. |5j to extract 
the infrared divergences. 

Our results have been obtained neglecting contributions of order aq/uMi and higher (q 
is the momentum transfer and M% is the mass of A). They cover both neutral and charged 
A and are reliable up to 0.5% or better in HSD. Furthermore, they provide a useful result for 
charm decay experiments with several thousands of events. For higher statistics experiments 
it will be necessary to incorporate aq/nMi contributions. 

II. VIRTUAL RADIATIVE CORRECTIONS 

In this section we shall discuss the virtual radiative corrections, up to order a and 
neglecting terms of order aq/irMi, to the DP of the hyperon semileptonic decay 

A^B + t + ut, (1) 

with A polarized. Our results will be especialized to exhibit the angular correlation s\ ■ 1. 

A convenient procedure to get such results is that of Ref. 0. So, in this paper we adopt 
the same approach, the same approximations and the same conventions of this reference. In 
this way, p\ = (Ei, pi), p 2 = (E 2 , p 2 ), I = (E, 1), and p u = (E v , p u ) will be the four-momenta 
of A, B, £, and i/£, respectively. M\, M 2 , and m will denote the masses of the first three 
particles. We shall assume throughout this paper that m v = 0. p 2 will denote a unit vector 
along P2, etc. We shall make our calculations in the center-of-mass frame of A. In this case, 
P2, I, and p v will denote the magnitudes of the corresponding three-momenta. There will be 
no confusion because the expressions obtained will not be manifestly covariant. Because we 
want our results to exhibit the correlation §1 • 1, it is convenient to choose the z-axis along 1 
and not along p 2 as in Ref. [§. 

At this point, it is convenient to mention that it is not possible to obtain the virtual RC 
of our present case by using the final virtual RC given by dTy of Eq. (15) of Ref. [[|. This 
is because in that equation the correlation §1 ■ p 2 was singled out after the integration over 
the azimuthal angle <j)£ of 1 was performed [in such Eq. (15) d(f>£ is still present in the phase 
space. As it appears and because of the choice of z-axis its integration amounts a 2tt factor]. 
Therefore, all the terms in dTy containing the product §1 • 1, and which appear before the 
integration over d<p£ is performed, have been transformed leaving the correlation §1 • p 2 only. 
There is no way to recover the §1 • 1 terms from this Eq. (15) of Ref. So, for obtaining 
our dTy, exhibiting the correlation Si ■ 1 only, we have to take a few steps back before that 
Eq. (15). 



3 



Our calculation now starts at the point where the scalar products • 1 and ■ p 2 appear 
for the first time, that is, at Eq. (13) of Ref. M, namely, 

^ \m v \ 2 = \y.\k\ 2 4eK 2 ■ (2) 



spins spins spins 



Here My is the sum of the order zero amplitude Mq, corrected by the model-dependent 
part of the virtual radiative corrections through the modified form factors f[ and g[, and the 
model-independent amplitude M v of such RC. Mq and M v are given explicitly in Appendix 
A [see Eqs. Q) and (|A3|)]. 

The two terms on the right-hand side of Eq. (|2]) arise after the spinor ua (pi) of the 
polarized hyperon A is replaced by £(si)-ua (pi) in the amplitude My. £(si) is the spin 
projection operator of A given in Eq. (4) of Ref. |4j]. 

With Eq. (H) we can express the differential decay rate dTy as 



, r dE 2 dE <m e dcj) 2 

dV v = — -g M 2 mm v 

(2tt) 



\H IKf-^Y, \ M v 

spins spins 



= dY'y-dY K y ) . (3) 

Notice the variables in the phase space of this equation, they correspond to our new choice 
of the z-axis along 1. 

In Eq. (|3p, dT'y corresponds to the first term within the square brackets. It can be 
identified with the differential decay rate with virtual radiative corrections of unpolarized 
A given in Eq. (10) of Ref. and, therefore, there is no need to recalculate it now. dTy 
corresponds to the second term of the square brackets of Eq. (H) and it contains the scalar 
products §x • 1, §i • p u , and s x ■ p 2 . The scalar product Si ■ p u can be expressed in terms of 
§! • 1 and §! • p 2 by three-momentum conservation. In this way, only Si ■ 1 and Si ■ p 2 will 
appear m dT { y>. Now, we require that §i • 1 be the only scalar product present in dTy . This 
can be accomplished by noting that the most general form of §i • p 2 depends on §i • 1, 1 • p 2 , 
and cos 02- The terms directly proportional to this cosine drop out after integration over 2 
from to 27r. This fact allows us to use the replacement 

§i ■ p 2 -> (§i • i)(l • p 2 ) = §i ■ iy , (4) 

in dTy ' , leaving us with an expression which only contains the correlation §i • 1. In Eq. (^) 
Uq is the cosine of the angle between p 2 and 1 when the real photon is not present. In the 
CM of A, Dq is given by 

V0 = 2p~l ' (5) 

where 

El = M 1 - E 2 - E, (6) 

is the neutrino energy (in this three-body decay). 

With these considerations in mind, we can express the dTy of Eq. (|3|) as 
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m G% dE 2 dE dflp dd)2~* r f oc , .. , ... ,,. 

y = -t (2tt) 5 1 r + * ( 10 + 10 } 



B'^-{B' 2 <t> + B' 2 \ 

7T 



This dr^ is the DP with virtual radiative corrections up to order a (and neglecting terms 
of order aq/irMi), leaving E 2 and E as the relevant variables and with only the angular 
correlation §i • 1 explicitly exhibited. The integration over 2 only amounts the factor 2ir. 
Now, A' Q , A'±, A" are given in Eqs. ( |Bl|) -(p3|) of Appendix B. The new terms Bq, B 2 , B' 2 ' are 

B'l = Q 6 Ep 2 y + Q r El, (8) 
B' 2 = -D 3 E°J + D,E (p 2 y + I) , (9) 
B'^ = D A E(p 2 y + l). (10) 

In these equations the coefficients Qq and Qj are long quadratic functions of the form 
factors. They are given in Eqs. (B6) and (B7), respectively, of Ref. 0. For completeness we 
repeat them too [see Eqs. ( [BPD and ( [BlUjJ ]. D 3 and D A depend on the leading form factors 
f[ and g[ and they are given in Eqs. ( |B13| ) and (P14j ), respectively. The primes on A' Q , A[, 
A", Bq, B' 2 , and B 2 indicate that these terms contain the model-dependence of the virtual 
radiative corrections through the leading form factors. The model-independent functions </> 
and 4>' are 

(j) (E) = 2 f — arctanh/3 — 1 ) In — — — (arctanh/5) 2 + — L 
\p Imp p 



+ iarctanh/?-il + | 7r2//5+ i ln(M2/m) (11) 
+ ^arctanh^ g + < ^ ln(Mi/m) nT ^„ < -I -I , 




(j)' (E) = (p-^J arctanh/5. (12) 

NDH corresponds to neutral decaying hyperons and CDH corresponds to charged decay- 
ing hyperons. The infrared divergence appears in the function cf) (E) through the parameter 
A. Let us mention that dYy of Eq. (0) contains two infrared divergent terms. The first one 
appears in A^cj) of the spin-independent part and the other one appears in B' 2 <f) of the spin- 
dependent part. Both terms will be canceled by their counterparts in the bremsstrahlung 
contribution. 

Let us close this section by comparing the dT v of Eq. (0) with the dT v of Eq. (15) of 
Ref. [Q. In spite of minor differences in their phase space factors, we can see that their 
spin-independent parts are the same. This is not the case for their spin-dependent parts. 
We can appreciate that the coefficients Aq, A' 2 , and A 2 , which appear in the spin- dependent 
part of dT v of Ref. J4|, have changed to the coefficients Bq, B' 2 , and B 2 , respectively, of dT v 
of Eq. (0). We observe in Eqs. (@)-(0) t na t in the B coefficients y always appears as a 



factor of P2, while for the A coefficients of Ref. Q yo always appears as a factor of I. This 
latter observation may induce us to think into the possibility of obtaining the dYy of Eq. ([/]) 
from the dYy of Ref. M by simply interchanging p 2 with I. Unfortunately this rule does not 
work because under such an interchange the A coefficients do not lead to the B coefficients 
and, thus, we can not obtain Eq. ((7|) directly from the final dYy of Ref. 
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III. BREMSSTRAHLUNG RADIATIVE CORRECTIONS 



In addition to the virtual RC the bremsstrahlung contributions must be calculated to get 
the complete RC to the DP of polarized decaying hyperons. In this section, we shall obtain 
them, in the same order of approximation as the virtual RC, both for the TBR and for the 
FBR. First, we shall define those regions and next, we shall proceed to the calculations. 
As is discussed in Appendix A, these corrections are model-independent by virtue of the 
theorem of Low ||. 



A. Kinematics, TBR, and FBR 

The DP in the variables E and E 2 is the kinematically allowed region of the four-body 

decay 

B + £ + v e + j, (13) 

where 7 represents a real photon with four-momentum k = (o>, k). The DP can be seen as 
the union of the TBR and FBR, each one defined presently. 

The TBR of the DP is the region where the three-body decay ([3]) and the four-body 
decay ( |T3"D overlap completely. The energies E and E 2 satisfy the bounds 

E 2 ~ <E 2 < E+, (14) 

and 

m<E<E m . (15) 
Here, E 2 (E 2 ) is the upper (lower) boundary of the TBR given by 

ei = \(M,-E±l) + nM M } E±lY (16) 

where E m is the maximum energy of the charged lepton 

Ml - Mf + m 2 
Em = 2M, • (17) 

The FBR of DP is the region where only the four-body decay ( |T3"D can take place. The 
energies E and E 2 in this region satisfy the bounds 

M 2 <E 2 < E 2 , (18) 
m < E < Eb, (19) 

where 

(Mi - M 2 f + m 2 
Eb ~ 2(M 1 -M 2 ) ' (20) 

Both regions TBR and FBR are depicted in Fig. 1 of Ref. || where more details about these 
regions can be found. 

We can now proceed to the calculation of the bremsstrahlung RC to the DP. First, we 
shall do this for the TBR and next for the FBR. The complete bremsstrahlung RC to the 
DP are obtained by simply adding the results of the TBR to those of the FBR. To obtain 
the complete RC of each region, we must also add the virtual RC of Eq. (0). 
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B. TBR bremsstrahlung RC 



Here we shall obtain the bremsstrahlung RC restricted to the TBR and with the angular 
correlation §1 • 1 explicitly shown. We shall follow the procedure introduced in Ref. [|J. 
Therefore, we start from Eq. (31) of this reference, 



El^ 



spins 



E 

spins 



\M' B 



spins 



M 



B 



(21) 



This equation is the square, summed over spins, of the bremsstrahlung transition am- 
plitude Mb of the four-body process (13) with the spinor ua(pi) replaced by Y*(si)ua(pi) ■ 
The explicit form of M B is given in Eq. ( |A5| ) of Appendix A. Equation (|21~1) enables us to 
express the bremsstrahlung differential decay rate as 



dT T B BR 



M 2 m m v d 3 p 2 d 3 l d 3 k d 3 p v 



(27T) 

dT' B - dT 



E 2 E 2uo E v 

W 

B ) 



E Wb\ 5\pi-p2-l-p„-k) 



spins 



(22) 



while dY 1 ^ contains 



where dT' B contains the first term of Eq. ( pTj) and is independent of §i 
the second term of this equation and is spin- dependent. dT' B is readily identified with the 
bremsstrahlung differential decay rate for unpolarized hyperons given in Eq. (33) of Ref. Q 



as 



dT' n = dT% + dT a R + dT 



B B ' B 

where, in accordance to Eq. (34) of this reference, 

a 
— < 

7T 



Bi 



dr% = -dn' {a[ [i (k/x) + c + d] + c 2 ) 



The phase space factor in Eq. (24) is now 



a a Gy dE 2 dE dQgd(j) 2 



-dVl' 

7T 



7T 2 



(27T) 



2M X , 



(23) 



(24) 



(25) 



instead of the dfl of Eq. (35) of Ref. ||]. The term I (k/X) fully contains the infrared 
divergence while C and C\ are the finite terms that come along with it. An alternative 
approach to extract the infrared divergence and the terms C and C% is the one of Ref. || , 
which was applied to the RC of DP of Kf 3 decays. This procedure was followed in Ref. 0] 

to extract the infrared divergence of the spin-dependent part dT^ B of its Eq. (32), and it 
can also be employed to extract the infrared divergence of the spin-independent part d T' B 
without any difficulty In Appendix C we give a brief review of this, here we only need the 
result of Eq. (56) of Ref. 0], which is 



I {E,E 2 )=I {k/X) + C + C 1 . 



(26) 



Here Iq (E, E 2 ) is the infrared- divergent integral given by Eq. ( |C6| ) of Appendix C. With 
this equation and the Eqs. (36)-(38) of Ref. H for C 2 , dY a Bl and dT b B , respectively, we can 
express the dT' B of Eq. (^), with some minor rearrangements, in the compact form 
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dT', 



B 



-eft)' A'Jo (E, E 2 ) + ^- dx dy d(f> k 
7r Air J-x J -i Jo 



with 



and 



\M 



„ a _ /3 2 (1 - rr 2 ) 
(1-Px) 2 



\M 



ll\2 



E„ 



ED (1 - px) 



Do 



D\E + D 2 lx 
D 



D 1 [u + E(l + Px) 



\M'\ 2 + \M"\ 2 



+ D 2 p u • 1 + k (E + u) - k 



m 



m 



E{l-(3x) 



E(l-/3x) 



(27) 



(28) 



(29) 



In these last three equations, x = 1 • k and y = 1 • p 2 are the cosines of the polar angles of k 
and P2, respectively whereas is the azimuthal angle of k. Furthermore, E v = E® — uj and 
D = E^ + (p 2 + 1) • k. The coefficients D\ and D 2 depend on the leading form factors and 
they are given in Eqs. ( pil|) and (pi2|) of Appendix B. Equation ( p7]) is also the explicit 
form of Eq. (55) of Ref. @. 



Once we have dT' B of Eq. (22), we can turn our attention to the spin-dependent part 
dr^g of this equation. In order to compute it we shall start from Eq. (43) of Ref. [[|, namely, 



Bi 



(30) 



where dY l B contains 



spins 



M, 



(s) 



+ 2Re 



M, 



and dT 1 ^ contains Zlspms 

dT l B contains the infrared- divergent terms as well as many infrared-convergent ones. dT l B 
is infrared-convergent only. To compute d T l B we follow the procedure of Ref. M to extract 
the infrared divergence (see Appendix C). According to this and using the explicit form of 



E 



spins 



given in Eq. (44) of Ref. 



] , we can write d T B 



as 



dT\ 



-dn'— lim 

7T 8lT A^O J\2 



d?k d^p u a . 
dtq— — — d [pi - p 2 



x 



u E u 

D 3 §i • IE% + DSi ■ P2E + ^§1 • IE 

l-e 



l-p v -k) 



+ D 3 s x ■ \u + DSi -kE}x^2 (j 1 



Pi ■ e 

k pi ■ k 



(31) 



The infrared divergence is contained in the first three terms within the square brackets, 



the remaining two terms are infrared-convergent. r\ = (p v + k) 2 is the invariant mass, which 
in the CM of A is given by 77 = 2p 2 l (y — y). In the TBR, ?] max = 2p 2 l (y + 1) and ?7 m i n = A 2 , 
with A 2 — > 0. The coefficients D% and D4 depend on the leading form factors and they are 



given in Eqs. (|B13| ) and ( |B14j ) of Appendix B. e M is the polarization four-vector of the real 



photon. 

In order to express d T B in terms of the correlation §1 • 1 we have to transform the scalar 
products §! • p 2 and §1 ■ k of Eq. (|3T|) in terms of §1 • 1. This can be achieved by using the 
substitutions of Ref. , 
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Si • p 



(§1 ■ i) (i ■ p) 



(p = P2, k, p„) 



(32) 



In Eq. (^Tp, with these substitutions, we can separate the infrared-divergent terms from 
the infrared-convergent ones as 

dT l B = -dQ% ■ l \b' — lim /"^ d v —^5 4 (pi- P 2-l-Pu- k) 
it [ 8tc a^o J\2 uj E v 



m 



M 2 



2pi-l 

Px-kl-k (I- kf (pi • kf 

d 3 kd 3 p v 4 
a?7 ^—o [pi - 



H lim , 

87T A->0 J A 2 



X 



zp 



P2 

p 2 1- 



l-p v -k) 



lk 



(33) 



(i-/3i-k) 

Here B' 2 is given in Eq. (||). In the first integral the sum over polarizations indicated 
in Eq. (31) was performed in covariant form, while in the second one it was performed by 
using the Coester representation M. 



The first integral in Eq. (|33"D can be identified with the divergent integral Iq of Eq. 
The second one can be put in the convenient form of Eq. (38) of Ref. by performing the 
integration over the 5 function and leaving y as the integration variable instead of 77. In this 
way, we can express d T B finally as 



dr R = -dn'si ■ i (b'i (e, e 2 ) + — [ dx f v ° d y r , 

71 I 4-7T J-l J-l JO 



D 



pi 



D d 1 



Ix 
D 



p 2 (1 - x 2 ) 
(l-(3xf 



(34) 



The other term of Eq. (|30|), dT l B , does not need any calculation. We can take the result 
of Eq. (57) of Ref. f|] and, with only minor changes in the phase space factor, we can adapt 
it to our case. After the application of rule (|3^) in such Eq. (57) is performed, we obtain 



dT B = -dQ si ■ 1— / — / dy / d(f) k — 

n 4:71 J-i 1 — Bx J-i Jo D 



x { D 3 E„ 



00 

-I - (E + u)x + 



m 



x 



+ Da 



uj + (1 + (3x) E - 



E 1 
2 1 



E 1-px 



px_ 

(p 2 y + I + u)x) 



(35) 



At this point we can collect our partial results to get our first main result, namely the DP 
of polarized decaying hyperons with radiative corrections up to order a, neglecting terms of 
order aq/nMi, and restricted to the TBR. It can be set as 



dT 



TBR 



TBR 



d Ty + d T B 
dT v + \dT' R - (dT R + dT^ 



(36) 
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where dYy is given in Eq. (|7]), <ir™ R is given in Eq. (|22|), and dT' B , dT l B and dT l B are given 
in Eqs. (p7f), (0), and (|35|) , respectively. The integrations over the three-momentum of the 
real photon in these last three equations are ready to be performed numerically (but they 
can be performed analytically also, as we shall see in the next section). 

The result of Eq. (^) can be compared with the corresponding result of Ref. ||] [the sum 
of Eqs. (15), (33), (51), and (57) of this reference]. We observe that the spin-independent 
parts dT' B are the same, although our present dT' B is expressed in a more compact form, 
while the spin- dependent parts are different. The results of Ref. 0] show that the correlation 
§i ■ p 2 is mixed with two others, §i • 1 and §i • k [see Eqs. (51) and (57) of this reference], 
while our present results exhibit only the correlation §i • 1. If we insisted in transforming 
the correlations s x • 1 and §i • k of Ref. @] as functions of Si ■ p 2 by using the substitutions 
of Eq. (D6) of this reference, the result would not give the same coefficient of §! • 1 as right 
here above. 



C. FBR bremsstrahlung RC and complete RC 



We shall now calculate the bremsstrahlung contribution of the FBR and afterwards 
we shall obtain the complete RC to the DP, with the addition of the TBR and virtual 
contributions. 

The calculation of bremsstrahlung in the FBR is relatively simple because the events in 
this region have the same amplitude M B of Eq. (|A5|) and it is infrared-convergent. We have 
to change only the upper limit of the integrals over the variable y of Eqs. (12 



©, and 

(j35|), which becomes one now, and to change the previously infrared-divergent integral Iq of 
these equations with I p , 



0, 



'OF 



OF 



hi 



' 3/0 + 1 ' 
,3/0-1. 



(37) 



with 



6*0^ = 4 (— arctanh/3 — 1 



(38) 



Iof is no longer infrared- divergent because in the FBR the photon has a minimum energy 
which is nonzero. It can be easily calculated from Eq. (|U4j) . The invariant mass rj of this 
equation must be now integrated from a minimum value rj m m = 2p 2 l (3/0 — 1) to a maximum 
value ?7 max = 2p 2 l (3/0 + 1). 

With these changes we can write the differential decay rate corresponding to the FBR 

as 



dT F B BR 



dT'i BR 



dY 



(s) FBR 



(39) 



with 



dT'i BR 



-dQ! 

7T 



A'J QF (E,E 2 ) + ^ 



dx / dy 



-1 



2n 



\M'\ 2 + \M'f 



(40) 



and 
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dr WFBR = a rffi / §i . J I b , j qf ,^ + P2i f 1 dx f 1 dy f 2w 

ix \ " An J-i J-i Jo 



IMT + 



M 



IV 



• (41) 



In Eq. (gg), |M'r and |M"f are given in Eqs. (|28| ) and (p9[), respectively. Equation ( f4"T|) 



is the sum of Eqs. ( |34]) and fl35|), after the above changes are performed. |M'"| and 
are 



\M 



W| 2 /3 2 (1-^ 2 ) 



(1 - 



D 3 l + D 4 Ex 1 



(42) 



M IV 


2 1 






DE 





m 



Da 00 



m 



1 



E l-px 



E 1-PxJ 
;i + /3x) e) {p 2 y + l+u) 



(43) 



Equation ( p5| ) is the FBR-contribution to the RC of the DP. It can be added to the 
TBR-contribution to obtain the complete RC to DP of polarized hyperons within the ap- 
proximations mentioned before. This completes our first main result of Eq. (^) by including 
the emission of all real photons allowed by energy-momentum conservation. It is displayed 

as 



dT = dT TBR + dT F B BR 



(44) 



with <ir TBR and dT™ R given in Eqs. (|36D and (0), respectively. The integrations over the 
photon variables are ready to be performed numerically. 

Let us close this section by mentioning that all the integrals which arise in the two 
regions of the DP can be performed analytically. Because of this, we can obtain completely 
analytical results for the RC of DP. We shall do this in the next section. 

IV. ANALYTICAL INTEGRATIONS 



In this section we shall perform analytically the photon three-momentum integrals con- 
tained in Eqs. fl2"T|), (|3"4"D, (f5j), and fl4"T| ) to obtain an analytical expression for the RC 
to the DP restricted to the TBR first and for the total DP later. 



A. TBR Analytical form 



The k-integrals corresponding to the TBR of the DP are those of Eqs. ( p7|) , (|34]) , and 
(|35|). They can be performed analytically by following the procedure of Sec. V of Ref. 0, 
where the RC of the DP of unpolarized hyperons were obtained. Fortunately much of the 
work has already been advanced. The integrals that concern us now can be expressed in 
terms of the functions 8i, i = 2, ... ,9 given by Eq. (99) of that reference, and in terms of 
the 8j, j — 10, . . . , 16 given by Eq. (46) of Ref. In this last reference the RC include all 
the terms of order aq/nMi, which are dropped here. We can express the analytical form of 
Eq. (|3) as 
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a 



dY' B = -dQ' A[I + {D 1 + D 2 ) 



71 



D 2 {9" + 9 IV 



(45) 



This equation is equivalent to dwB of Eq. (92) of Ref. 0. Similarly, the expressions for 
dY l B and dY l B of Eqs. (0) and (|35| ) become 



dT\ 



r , - -dQ' §x ■ 1 \b' 2 I + D 3 pi + D 4 p £ 2 

7T 1 



(46) 



dT l l = -dtt'sfl D 3 pi + D 4 pl 



71 



(47) 



In Eq. fl45|) the 0j functions are contained in the functions 9', 6", 9"' and IV given by 
Eqs. (85), (86), (90), and (91) of Ref. ||, respectively. The functions p\, i = 1,...,4 in 
Eqs. (|46|)-(|47|) can be expressed as 



A 



P3 



Pi 



2 

P2E 2 



E 



1 2 + 20 3 - 4 



1 02" 



0o + 



El 

2 



E(E + E°)(l-f)9 2 - (3-P 2 )^ + EE° 



-3)03-204 
£_ 2 

T 

,2 



2 

+ - [2E-E")9 



P2 

2 
+ 



— 8 + ^09 — 



2£ 

?9 



10 



'15 



m 



2-/3 2 + ^ 
£ 2 



02 + 



-^m 2 + p 2 ^o 



^-^) — -EEl- P2 ly, 



-(3E 2 + 2Z£° 
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—0 6 + - (3E + p 2 (3y ) 7 - E9 8 + -0 9 - -/ 2 

P2 l (l - /? 2 ) n + 2 P2 /0 12 - P2 /0 13 - ^014 - J015 - t^0 



(48) 
(49) 



(50) 



(51) 



With Eqs. ©-© all the integrals over k in Eqs. fl27|) , (34]), and (^) have been 
expressed in an analytical form. We can obtain now the bremsstrahlung differential decay 
rate o?r™ R of decaying polarized hyperons with the photon integrals expressed analytically. 

) the analytical forms of dT' B Eq. (45) and of dT B \ which is the sum 



Substituting in Eq. ( 
of Eqs. ( f4"6] ) and (|4"7|), we obtain the analytical form of <fT™ R 



dT T B BR 



a 



dQ' {A[I + p! + D 2 ) (0' + 9"') + D 2 (9" + 
- h ■ 1 [£? 2 / + D 3 (Vi + P3) + ^4 (p £ 2 + Pi)} } • 



iv 



(52) 
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We are now in a position to obtain our second main result in this paper: the analytical 
RC to the DP of decaying polarized hyperons to order a and neglecting terms of order 
aq/irMi. This result comes from the addition of the virtual RC, dTy of Eq. (0), and of 
^pTBR Q £ (|52|) , It can be put compactly as 



dY 



TBR 



Gl dEo dE dVto 



a 



2M 1 A' +-$!-§!■! 



7T 



2 (2nY 

Here A' and B' ' are the same as Eqs. ( pi|) and 
$x = A\ (0 + J ) + Aftf + (D x + D 2 ) (9' 



B" + -$ 



7T 



(53) 



respectively. $i and $2 are 

"') + d 2 (e" + # IV ) , (54) 



= B' 2 (0 + Jo) + S^' + ^3 [P\ +Pl)+ D A U 2 + pi 



(55) 



The coefficients A' x , A", B' 2 , and 5 2 ' are given in Eqs. (§2|), (|B| 
tively. Di, i — 1, ... ,4 are given in Eqs. ( Bll|) -( |BT4l) . The functions 



and (|10"1), respec- 
, and Iq appear in 



were given in Eqs. 
and (94) of Ref. [| 
term has to be ^9 5 



Eqs. ([TT|), (0), and QCT1|), respectively. The new model-independent functions pf, i = 1, . . . , 4 
H)-((51|). The sums & + 9"' and 0" + 9 W appear explicitly in Eqs. (93) 
respectively. We have corrected a misprint in that Eq. (93). Its fourth 
rather than — \9^. For completeness, we reproduce these two sums in 
Appendix B [see Eqs. (|B39p and (|B40| )1. 

Because the infrared divergence, which appears in the virtual part dTy through the 
function 0, cancels out with its bremsstrahlung counterpart, which appears in Jo, the sum 
4> + Iq is no longer infrared-divergent and, therefore, o?r TBR of Eq. (|SB|) is infrared-convergent. 

^pTBR - g ^ e result corresponding to dY of Eq. (101) in Ref. f|]. Comparing both results, 
we can observe that the spin-independent parts are the same, but the spin-dependent parts 
show important differences. In Eq. (101) of Ref. f| we have, within the square brackets that 
accompany the correlation §1 • j>2, the terms A' Q ' and $2, while in Eq. (|53|) the corresponding 
terms accompanying §1 ■ 1 are Bq and $3. They are different. In fact, we may notice that in 
$2 °f Eq. (^) only the ^-functions, i — 2, . . . , 16 appear, while in <E> 2 of Eq. (103) of Ref. |4| 
also the //-functions appear [see Eqs. (91)-(95) of this reference]. Because of this, the RC to 
§! ■ 1 and §! ■ p 2 correlations are quite different. 



B. FBR analytical form 



The k-integrals of the FBR are those contained in Eqs. fl40| ) and fl4"T|) and they can be 
performed analytically, too. Because <ir'J rBR of Eq. (flOP has the same form as the corre- 
sponding dY' B of Eq. ( |27| ) of the TBR, we can follow the same procedure of the Sec. V of 
Ref. to calculate the analytical form of this dT'^ BR . The result has the same structure 
as dT' B of Eq. 



dT'i BR 



-dtt' 

7T 



A[I 0F + (D 1 + D 2 



+ D 2 



9" +9™ 



(56) 



with 
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9' F + e'" 



-El (i - p 2 ) e 2F + (eI - 1 -±^e 



P2I 
2 



2 — 2 — 




™-E$ 1 1 " 



and 



2 



O f - + El + /3p 2 y Q ) # 3 f + (El + E) 6 AF 



w 



5F 



(57) 



(58) 



These last two equations also have the same structure as 9' + 9"' and 9" + # IV , which 
appear in the analytical form of dT' B of Eq. fl45|) . The difference between dT' B and dT' B 
of the TBR lies in the integrals Iq, Iq F and in the functions 0j. In the FBR these functions 
change to 9i F because the variable y must be integrated from —1 to 1 and not from — 1 
to yo as in the TBR. The new set {9 if} ,i = 2, . . . , 16 is explicitly given in Appendix B 
and 9of is given in Eq. (|38"D . We can compare with the 9f of Ref. where the RC of 
DP for the FBR were calculated up to order aq/irMi. We cannot take readily the result 
of this reference because, according to our approximations, we would have to neglect all 
the terms of order aq/nMi in that result to obtain ours. We find the procedure of Ref. |§ 
more adequate and straightforward for our purposes. However, in order to check our results 
we have reproduced the Table I of Ref. ||. Our numerical evaluations coincide very well 
within the approximation of neglecting terms of order aqji\M\\ we shall not display here 
this numerical evaluation. 

In a similar way, we can see that the spin-dependent part of <ir™ R has the same structure 
as the corresponding spin-dependent part of (£T™ R . Thus, we get 

a 



dT 



(s) FBR 
B 



71 



-dtt' §1 • 1 



B' 2 I 0F + D 3 (p* F + p£ F ) + D 4 (p 2 £ F + pi F 



(59) 



with the functions pf F , i - 
(|5"T|). Explicitly, they are 

P2I 2 



4 having the same structure as the previous p\ of Eqs. (^8| 



Pif 



P2F 



2 

p 2 E 2 



— 1 ) 02F + 2#3_f — 9. 



'AF 



~E 9 ° F + 



1)9 



2F 



(f3 2 - 3) 9 3F 



2,9 4 F — 



75F 



(60) 
(61) 



P3F = ^\e(e + El) (l - p 2 ) 9 2F - (3 - P 2 ) ^ + EEl 



+ \e 2 (l + /? 2 ) 9 AF - l -(E + 2El) 9 5F - ^ 
+ i(2E-E£) 9 7F -\(E-E[ 



2E 

9&F + ~^@9F 



>3F 



-I 9 WF — -@15F 



(62) 



PiF 





r e°~ 








jm 2 


2-/5 2 + f 


92F + 


- 7 -m 2 + p 2 lyv 


9 3 F 



(3 - (3 2 ) ®- - EEl ~ P2lyc 



'AF 



-f3E 2 
2 H 



21E°, 



'5F 
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- + - (SE + P2 (3y ) 6 7F - E8 8F + -8 9F - -l 2 6 10F 

- p 2 l (l - /3 2 ) 6 nF + 2p 2 Wi2F - P2W13F ~ ^0 14 F - ^i5F - ^JjAeF j 



(63) 



From Eqs. ( |56| ) and fl5"9|) we obtain the analytical bremsstrahlung differential decay rate 
^-pFBR Q f decaying polarized hyperons corresponding to Eq. (j39j), 



with 



dT F B BR = -dD! 

7T 



$1F - §1 • 1 $ 



2F 



(64) 



$ 1F = A;/ 0F + (A + D 2 ) (^ + 0£') +D 2 (9 F + 6™ 



(65) 



(66) 



$4 = ^/ 0F + D 3 (p/ F + p/ F ) + £> 4 (p 2 £ F + p/ F ) . 

At this point we complete our second main result. The addition of dT™ R of Eq. 
and <ir TBR of Eq. (|53| ) gives us the complete analytical RC to the DP of decaying polarized 
hyperons to order a and neglecting terms of order aq/irMi. This complete result can be 
expressed compactly as 



dT 



TOT 



G$ dE 2 dE dVL t 

"2 (27T) 4 



a 



2Mx 4, + - ($ ! + $ 1F ) - s x • 1 

7T 



(67) 



Here A' Q , B'J, $ x , $i F , and $^ are given in Eqs. (jm]), (§), fl54|), (H), and ©, 

respectively. 

From Eq. ( |6"TD we can obtain easily Eq. (|53| ) to the RC of DP with the TBR only by 
dropping $ 1F and <& 2F - It is this Eq. ( |67D which must be used to obtain, in HSD, totally 
integrated observables, such as the spin-asymmetry coefficient of the charged lepton. We 
shall calculate this asymmetry coefficient in the next section, allowing for the possibility that 
real photon emission be discriminated experimentally via energy-momentum conservation. 



V. SPIN ASYMMETRY COEFFICIENT a t 

In this section we shall obtain the RC to the spin-asymmetry coefficient of the charged 
lepton ct£. We shall consider the two cases discussed all along, namely, that bremsstrahlung 
photons not be discriminated at all or that directly or indirectly the photons belonging to 
the FBR be eliminated from the experimental analysis. We will discuss the former case first 
and afterwards we will discuss the latter case. As we shall see in the numerical evaluation 
of the next section, an appreciable difference can be observed between these two cases. 

ct£ can be calculated from the total DP of Eq. (|67l). This equation can be used to get 
the quantities which appear in the definition of «£, 

N + - N~ 

a e = 2— — . 68 

N+ + N- v ' 
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Here iV + (N~) denotes the number of the emitted charged leptons with momenta in the 
forward (backward) hemisphere with respect to the polarization of the decaying hyperon. 
With those numbers calculated, we may express as 



B 2 e +(a/7i) (o/ + 



l 2F 



a J — ^ 7 ' ; T7~ —T- (69) 

Bi + (a/7r)(ai + aur) 



Here 



E m rE^ 



B l 2 = I "' /J B%dE 2 dE, (70) 
$/d£ 2 d£, (71) 



r E B f E 2 



4f = / / H F dE 2 dE, (72) 

J Ah 



2 



fii = / A(,d£ 2 dE, (73) 



a x = § 1 dE 2 dE, (74) 



E ~2 

rE B rE~ 

a 1F = / <S> 1F dE 2 dE. (75) 

In these integrals, B%, ® e 2 , ¥ 2F , A' , $i, and $ 1F are given in Eqs. ®, (H), ©, TO, (H), 
and flBUD, respectively 

In Eq. fl6"9|) we have attached an upper index T to denote that the asymmetry coefficient 
includes the total DP of the real photons. The contributions of the TBR to the RC of ae are 
given by the terms a 2 and a x , while the contributions of the FBR are given by the terms a 2F 
and a\p. We can now rewrite aJ to comply with our approximations, i.e., in such a way that 
only the terms of order a, neglecting terms of order aq/nMi, appear. The corresponding 
expression is 



a e = a 



1 + - 



a I a 2 + a 2F ai + dip 



vr V B' 2 (0) B x (0) 



(76) 



where a e Q is the spin-asymmetry coefficient of the charged lepton without RC. It is obtained 
from Eq. (K^) by dropping the terms proportional to a, namely, 

< = • (77) 

B 2 (0) and B\ (0) in the denominators of Eq. ( |76[) are the zero-order q/M\ approximations 
of the B 2 of Eq. (|70|) and of the B x of Eq. (|73|), respectively. Explicitly they are, 
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B l 2 (0) 



B 1 (0) 



B' dE 2 dE, 



A[ dE 2 dE, 



(78) 
(79) 



with B' 2 and A[ given in Eqs. (||) and ( |B2|) , respectively. 

The coefficient for at when only the TBR of the DP is allowed can be easily obtained 
now. All that has to be done is to drop a 2F and a\p from Eq. ( j76|) so that 



a, 



a I a 2 ai 



(80) 



We attached an upper index R to denote that the bremsstrahlung correction is restricted to 
the TBR. 

In Ref . [|J we only calculated the emitted baryon asymmetry-coefficient a b corresponding 
to af. In this reference it was assumed that the FBR photons were always discriminated. 
The contribution of these photons should be calculated and added to the results of this 
reference in order to get ana^ corresponding to the above aj. 

In the next section we shall display numerical evaluations that will allow us to compare 
our results with others available in the literature and, also, to appreciate the relevance of 
discriminating or not FBR photons. 



VI. NUMERICAL RESULT 

In order to compare the coefficients af of the TBR of the DP and aj of the total DP, 
we shall make numerical evaluations of them for several decays. These results will enable 
us to establish the relevance of the difference between af and aj in the study of HSD. We 
shall also compare them with other results reported in the literature. In Table I we give 
the values of the form factors used in the numerical evaluation of the coefficients af and 
aj for the decays n — > peV, A — > peV, S~ — > neu, £~ — > AeV, S + — > Ae + v, H~ — > Aez7, 
2~ — > S°ez/, S° — > E + e77, and A+ — > Ae + u. For this last decay we take the form factors 
of Ref. DTof. The sign of the form factor gi must be changed when the charged lepton is 



positive [|TTJ. In the radiatively uncorrected amplitudes the g 2 -dependence of the form factor 
was neglected along with the contributions arising from / 3 , g 2 , and as was done in other 
calculations in the literature. 

To evaluate af we use Eq. (|80|) and for aj we use Eq. (|76|). These equations involve 
the double integration over the energies E and E 2 . At this point is convenient to mention a 
technical aspect that we have to deal with in calculating the integrals 



i 

Je. 



In (y + 1) dE 2 , (81) 
and 

rE+ /,/„-!- 1\ 

(82) 
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which are contained in the infrared- divergent integral Iq of Eq. ( |C6|) and in the integral Iqf 
of Eq. (j37[) , respectively. The integrals of Eqs. fl8lD and (|82D diverge when ?/ — ► ±1, that is 
when E 2 = Ef. However, it turns out that the result is finite as was shown in Ref. [Q. To 
calculate it we can either follow this approach and implement it in the program to evaluate 
the a^s or we may neglect the points where yo — > ±1 . This last is equivalent to leaving 
out the boundaries of the TBR and the FBR of the DP. The numerical difference between 
these two alternatives is negligible. Here we follow the second alternative. 

In Table II we display our numerical results for the radiative corrections of the asymmetry 
coefficients af and aj . We compute them by taking the percentage differences (that is, we 
multiply by 100) 

r R,T r,t e / q o \ 

oa e = a e — a , {SS) 

where a l is the uncorrected spin-asymmetry coefficient of the charged lepton, Eq. (|77|) . 

In the second column of Table II we display the <5a R corresponding to the TBR of the 
DP, in the third column the Sa T corresponding to the complete DP are given, in the fourth 
column we give the results for 5a obtained from Eq. (23) and Table I of page 58 of Ref. JTH 
and, finally, in the fifth column we give the two values reported in Ref. |J. 

From this Table II we see that there is a very good agreement between our 5a T and 
the Sa of Refs. |H| and || . In both these references the FBR was included. The inclusion 



or exclusion of the FBR is appreciable, as can be seen by comparing the second and third 
columns, except for the decays n — > peV and £~ — > neV. In several instances the inclusion of 
the FBR contributions reduces the total radiative corrections, even to the point of making 
them negligibly small. It may be that the values in the second column are one order of 
magnitude larger than the corresponding ones in the third column. Therefore, in general, 
there is an important difference between af and aj. Experiments of HSD where no provision 
to detect the photon has been made must use aj to study such decays. But, if the photon 
is discriminated eliminating the FBR of the DP then the adequate coefficient to study HSD 
is af. 



VII. CONCLUSIONS 

In this paper we have obtained the radiative corrections to order a to the Dalitz plot 
of the semileptonic decays of polarized spin-1/2 baryons, neglecting terms of order aq/irMi 
and higher. Our main result, which is compactly given by 

dT = dT TBR + dT F B BR , (84) 

is presented in two forms: one in which the triple k-integration is ready to be performed 
numerically, Eqs. fl36f) and (@), and another one in which such integration has been an- 
alytically performed, Eqs. (|53| ) and (p7|). Since real photons may be discriminated either 
directly (by detection) or indirectly (by energy-momentum conservation) we have split the 
above two forms to cover this possibility. In case this photon discrimination takes place, 
Eq. ( |34D becomes 



dT TBR = dT + 



dT' B ~ (dT l B + dT 1 ^)] , (85) 
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where dTy is given in Eq. (0) and it corresponds to the virtual RC. dT' B , dT l B , and dT l B 
correspond to the bremsstrahlung RC and are given by Eqs. (p7|), fl34|) , and (^), respectively, 
with the photon integrals remaining to be performed numerically. The infrared divergence 
and the finite terms that accompany it have been explicitly and analytically extracted, 
however. They appear in the Iq of dT' B and of dT l B . The analytical result for Eq. (^5|) is 
given in Eq. (|53]). 

When real photons are not discriminated what-so-ever, then the DP with the union of 
the TBR and the FBR must be used, as in Eq. fl84l) . <fT™ R is given in Eq. (^), where the 
k-integration remains to be evaluated numerically, and the analytical result for Eq. flSj|) is 
given in Eq. (|6T|) . 

An important integrated observable is the charged-lepton spin-asymmetry coefficient ag. 
Using the analytical form we obtained the radiative corrections, through Eq. (|6T|), to this 
observable. The integrations over E and E 2 were performed numerically and the results are 
displayed in Table II. A systematic behavior of the RC to ae is observed. The contribution 
of the FBR bremsstrahlung may be as important as the RC from the TBR and even of 
opposite sign, in such a way that when no photon discrimination takes place the complete 
RC to oti may become almost negligible. In this Table we also compare with results reported 
in other references. This comparison is satisfactory. 

To our knowledge Eqs. (|53|) and (^) for the RC to the DP allowing for the TBR only 
and including the FBR are the only analytical expressions available in the literature. In 
Ref. || the relevant variables of the DP are also E and E 2 , but there only numerical results 
are given, which are compromised to particular values of the form factors. Our numerical 
evaluations coincide very well with the numbers reported there for the percent RC of the 
asymmetry coefficient aj for the decays £~ — > neV and A — > peV. We have also evaluated 
the RC for the process A+ — > Ae + u. The numbers obtained are presented in Table II and 
in this decay they are a good first approximation, which is already useful given that the 
experimental statistics are not yet too high. 

Our results are model-independent and are not compromised to any particular value 
of the form factors. All the model dependence of radiative corrections has been absorbed 
into the f\ and g\ form factors in our approximation of neglecting contributions of order 
aq/irMi. This is indicated by putting a prime on them. For hyperons our results are reliable 
up to a precision of around 0.5%. This precision is useful for experiments involving several 
thousands of events. For high statistics experiments involving several hundreds of thousands 
of events or for decays involving charm as A+ — > Ae + u or even heavier quarks our equations 
provide a good first approximation. To improve the precision of our formulas it becomes 
necessary to include aq/nMi contributions. Our results are valid both for neutral or charged 
decaying hyperons and whether the emitted positively or negatively charged lepton is either 
electron-type or muon-type. To conclude let us remark that in a Monte Carlo analysis the 
advantage of the analytical form is that the triple k-integration does not have to be repeated 
every time the values of fi and gi, or of E and E 2 , are changed. This leads to a considerable 
saving of computer time. 
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APPENDIX A: 



In this Appendix we give for completeness the amplitudes for the RC of the decay ([!]). All 
of them are also given in Ref. ||. The uncorrected transition amplitude Mq for process ([!]) 
is 



G 

M = —= [u B (p 2 ) Wft (pi,p 2 ) u A (pi)] [ut(l)OpV v (p„)] 



where 



V2 



(Al) 



M 1 



h{q 2 



Mi 



+ 



92(q 2 



Ouvqv + 



9s(q 2 

M 1 



75- 



(A2) 



Here M = 7 M (1 + 75) and q is the four-momentum transfer. 

The model-independent part of the virtual radiative corrections has the amplitude 



a 



M v = — [M <t>(E)+M Pl 4>' (E)}, 
where 4>(E) and 0' (E) were given in Eqs. (|TI|) and fll2"|) , respectively. M P1 is 



M 



in 



\mMj 



E \ G 

-j= {ubWxu a ] [ue^Oxv^ 



(A3) 



(A4) 



The model-dependent part of the virtual radiative corrections is absorbed into Mq 
through the definition of effective form factors f[ and g[. This fact is denoted by putting a 
prime on Mq. The bremsstrahlung transition amplitude Mb is obtained following the Low 
theorem M, 



M b = —^=- [u b W x ua] [uiOxv^ 
= M a + M b . 



I 



Pi ■ e 



I ■ k pi ■ k 



&Gv r _ TTT 1 

+ -j=- [UbWxUa] 



21 ■ k 



(A5) 



Within our approximation the Low theorem guarantees that no model-dependence appears 
here. 
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APPENDIX B: 



In order to make this paper self-contained, we reproduce here all the coefficients which 
appear in our final results. For Eq. (J7|) they come from Ref. [|J, 

A' = Q x EEl - Q 2 Ep 2 (p 2 + ly ) - Q 3 l (p 2 y + I) + Q 4 E° u p 2 ly - Q 5 p 2 2 ly (p 2 + ly ) , (Bl) 



A[ = D 1 EE° U - D 2 l (p 2 y + I) 



(B2) 



A'[ = DiEE®. 



(B3) 



The coefficients Qi, i = 1, . . . , 7 are given in Eqs. (B1)-(B7) of Ref. ||, respectively. For 
completeness we reproduce them here. 



Qi = Fi 



Q 6 



2E 2 - M 2 



x 



Mi 
E 2 

Mi 



1 

2' 



M 2 + E 2 



+ F 2 Fz 



Mi 
M 2 + E 2 



+ FiF 2 



M 2 + E 2 



Mi 



1 + 



M2 
Mi 



Mi 



E 2 
Mi 



2E 2 + M 2 



Mi 



G\ 



M 2 -E 2 



2 

G 2 G% 



Mi 
M 2 -E 2 



M 1 



+ G X G 2 
1 



M 2 -E 2 



E 2 
Mi 



Mi 



+ C1G3 



M 
M x 



_E 2 

Mi 



Mi - E 2 



Mi 



2 Ml}' 



F"l G\ F\F 2 G\G 2 FiF s 



Mi Mi Mi Mi Mi 



Mi. 



+ 



E 2 F 3 



Mi 



M 2 + E 2 



+ 



G1G3 



Mi 

Qz = Qi- IFl 



M, 



- 1 



G 2 G% 



Mi 
E 2 -M 2 



Mi 



Mi 
-2G? 



M 2 -E 2 



Mi 

E 2 + M 2 
Mi 



+ 24^ + M1Q5 



Mi 



M^Qr 



£/2 

Mi 



Mi 
Mi - £ 2 

M 

n 2 „2 ' 



g 

Mi 2 /' 



Q 4 = Q2 - 4 



F1G1 
Mi 



if_ \M 2 + £ 2 

Mf 



M 1 

E 2 - M 2 p 2 f3y 



g| tm 2 - £ 2 

Mf 



Mi 



- 2-t-J- + 2- 



Ml Ml 



M 1 

{( n G 2 - FiF 2 ) 
M 2 



Mi 

M x 



£ 2 + M 2 p 2 /% 



M x 



Mi 



+ F 2 G 2 



+ 2FiGi 



1 + ( 1 ^ 2 )£+t+ 



E 2 - p 2 py 



Mi 



+ F X G 2 



-1 + 

Mi 
E , p 2 /%o 



+ (l + /5 2 ) 



Mi M x 



Mi Mi 
2 / 



M x Mi Mi 
Gi_F 2 



Mi 



— F3G3 



m" l E 2 / 2 \ i£ P20Vo 

m|I v 1 ~m7~v 1 "^Jm7 + "mT 



(B4) 



(B5) 
(B6) 
(B7) 
(B8) 
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m 2 ( M 2 E 

+ — + — 

M X E \ Mi Mi 



F3G1 



m 2 / M 2 E 

-1 + — 

M X E \ Mi Mi 



(F 2 G 3 + F 3 G 2 



m? (Mx-E 2 -E 



M\E V Mi 



(B9) 



(Mi + M 2 ) {E 2 - M 2 ) 



+ G\ 



(M x - M 2 ) (£ 2 + M 2 ) 



+ 2FiGi 



MiE 

Mi (-Mi + E 2 + 2E) - m 2 



x 



Mi -2E-E 2 
E 

E 2 + M 2 
Mi 



MiE 

- GiF 2 

f m 2 N 
V M£, 



M X E 
+ FiG 2 



ff 2 -M 2 
Mi 

£2 + M 2 \ /Mi - 2E- E 2 



Mi / V 



G3F1 



E 2 -M 2 
Mi 



MiE 



The coefficients D~, j = I, . . . ,4 are 



D 3 = 2(f[g[ - g[ 2 ) 



D, = 2(f[g' 1 +g[' 



(BIO) 

(Bll) 
(B12) 
(B13) 
(B14) 



and ?/o were defined in Eqs. (§) and (|6]), respectively. 

The functions 9i which appear in Eqs. (^8|)- (|5T|) corresponding to the TBR of the DP 
are given by 



0; 



~{T l + + T i r 

P2 V 



(B15) 



where i — 2, ... ,16, and 



1 2 



In 



lT/3 



± 



(lT/3)(l + /3a±) 



1 — (3xq 
ln(l ±a ± ) 



± 



(l±x )ln(l±x ) 
(l±/3)(l-/?x ) 



(l±/3)(l + /3a±) 

1 ± a 1 * 1 , . ± , (xq + a 1 * 1 ) ln(±xo ± a 



(l + /3a±)(l-/?x ) 



(B16) 



1 



3 3 2(3 

1 - px 



{' 


r 1-/3 ■ 


- L 


1 — /3x 


- L 


"1 + (3a-~ 




1 - /3x _ 




. 1 + P . 




1 — (3xq 



+ L 



1+pa 
1 + P 



+ L 



l + {3a^ 



1 + 



In 



1 ~ Px 
1-/3 



In 



1 + /3a 4 

T+p 



(B17) 
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(x ± 1) ln(l ± x ) ± (1 ± a ± ) ln(l ± a ± ) - (x + a ± ) ln(±x ± a ± ), 



-i {(1 - z 2 ) ln(l ± x ) + (x T l)o ± + 1 - (1 - a ±2 ) ln(l ± a ± ) 
+ (x2-a ±2 )ln[±(x + a ± )]}, 



-/ + p 2 ± 



/3E°(x + a J 
1 + /Sa^ 



E0 _PE°(x + a*) 



(l+(3a^f 



(1 + Pa*) 



(l + /3a=F)5 



P2 - Z =F 



± 



± J 2 



f3E^(x + a ± ) 



l + (3a 



± 



Ji 



El (x + 



J 2 i 



l + f3a± 

-2(1 -p 2 + E° v x ) T El(x + a ± )/ 2 ± - E° v (x + a ± ) J±, 



-^(l-p 2 + E° u x ) + 



3(1 -p 2 ) | | pc± 



A(3l 



4/ 2 



(^) 2 K + « ± ) 2 r± 
4/ 2 (l +/3a ± ) 



±^ 3 



(£°) 2 (x + a 
4/ 2 (l + /3a±) 



±\2 



±x -3 



3£° + 3E°(Eg + Zx ) ± ^ 



4/5/ 



4/2 



Ji ± G ± Jt 



^(xl =F 1) ln(l Ti ) + | [(a T ) 3 =F l] ln(l Ta T )-| [4 + (^ T ) 3 ] h [t^o + a T ) 
+ i(l-^)(a^±l)-i( a ;o±l)[l-(a^) s 



T n = ^ [(1 - /fc ) J 4 " Ji] " (PE° V + / - p 2 )h + (l~ P2)h} , 

E° v (l - /3x ) J, + 2E° u x + 2(1 - p 2 ) - ((3E° + I - p 2 )h 
l'v, = -—El(l-xl), 



12 



1 

1 



2p 2 



El 



l + x 2 + 2a ± (x =F 1) ± a ± (x + a ± )(/ 2 ± ± ^ 



3E° U [2p 2 (l + y ) + 1(1 - x 2 )] - (E° v ) 2 (x + a ± ) 2 ( J 3 ± ± J 3 ± ) 
-2lEl(x + a ± )a ± (J 2 ± ±lf) ) 
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2/? 2 



2(l-p 2 + E „x )+pE°„(l-xl) 



+ 



3(l-p 2 + (3El 



2p 2 



+ 



+ 



■A 



2/3 2 2/ 2 
The following definitions are used in the above expressions 



x 



P2Vo + 1 
El ' 



I 



Ix = — arctanh/3, 

P 



If = In 



a ± + 1 



a± - 1 



a ±2 - 1 



J. 



In 



1 - (3x 



+ In 



1-/3 2 ' 
1-/3 



+ In 

a* + 1 



1 - (3x Q 



J 



-2 



a ± + x 
1 



Ja = ~; 



a ±2 — 1 a 1 * 1 + Xq 
1 1 



1 - 13 2 l-/3xo 



/3(E°) 2 (xq + a^) 2 a* (a ±2 - 1) 
~F .,o „ "F 



4/ 2 (l + /3a±) 2 ' 4(l+/3a±) 
The sums 9' + 0'" and 9" + # IV which appear in Eq. are 

^ (l " /^ 2 ) ^2 + U - 

2E - E?. „ 1 



0' + 9" 



P2I 
2 



2 A 2 



+ -9 5 H — / 

2 2 



2£ 



— 1/0 — — ~9q 

7 2 8 4E 9 



(B30) 



(B31) 
(B32) 
(B33) 
(B34) 
(B35) 
(B36) 
(B37) 
(B38) 



(B39) 



9" + 9 lv 



P2I r 
2 



9 -(E + E° u + Pp 2 y ) 9 3 + (E° u + E)9 A + W* 



(B40) 



The explicit form of the photon integrals corresponding to the FBR of the DP is 



>2F 



PP2 



b+ 



E' 



J 2 +-/ 2 -+/31n 



+ 



2h 

Eb-b+ 1 



(B41) 
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?3F 



6qf 
9 7F 

08F 



>10F 



7 11F 



7 12F 



'UF 



h 



In 



Q§f — - — 



P2 

1 

P2 
1 
2P2 



b + 


1 


H 


b- 


+ PP2 





l 1z£\ +l 1+£)-l 



b+ 



b+ 



a + I 9 + — a I 9 + In 



(l-a +2 )/ 2 + -(l-a- 2 )/ 2 - + ^ 



2^(/ 2 -+/3/i)-2^(/ 2 + + /3/ 1 )+2 



2 + p\ p 



b+ 



2/ 1 + ^(/?/ 1 + / 2 -)-^(/?J 1 + V 



6 9F = 24E + 2 



6(E° U -E)+ P(G F + G+)} h + 2(G F I 2 + G + F P 2 



+ 2P2 



b- 3 b+ h 



i 

2(/ 4 - h) 

PP2 

2{h - 2) 



2(a' 2 - a+ 2 ) - a" 3 /^ + a +3 I 2 + + In 



#13F — 0, 



PP2 



(2-a-I 2 -)(y )-(2-a + I 2 + )(y+) 



9 15F = 24E?, + 4l 



a y I 2 - a + y+I 2 + + 2p 2 (y^flf - (% + ) V 



6 mF = 2AE 2 (h - 2) + 8[(E^) 2 - 2E 2 (5 2 ]I 



+ 4:lp 2 



{y ° )2 iPh + I 2 -)- { 4^(Ph + I 2 



b- v x z ' ¥ 
where a ± , Ji, I 2 , I3 , J 4 are given in Eqs. ( [B31| )-( |B"3~3| ) and 

b ± = l + f3a ± , y^=y Q ±a ± , 



GP = T P (^Ea^+p^) ff. 



(B42) 

(B43) 
(B44) 

(B45) 

(B46) 
(B47) 

(B48) 

(B49) 

(B50) 

(B51) 
(B52) 
(B53) 
(B54) 

(B55) 
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APPENDIX C: 



In this appendix we give a brief review of the procedure of Ref. || to extract the infrared 
divergence and the finite terms that come along with it. This procedure can be adapted to 
our case by expressing the differential decay rate of the bremsstrahlung radiative corrections 
(1Tb in terms of the invariant mass i] = (jp v + k) 2 as follows 

d 3 k d 3 p v 



dT 



M 2 mm u dE dE 2 dVt^ 
(2tt) 8 4 



drj- 



u E v 



Y J \MBY5 A { Pl -p 2 -l- Pu -k). (CI) 



spins 



For the TBR and in the CM-frame of the decaying particle, r/ = 2p 2 l(yo — y), and 
Vmax = 2p 2 / (y + 1) and r/ min = A 2 , with A — > 0. Here A is a small mass assigned to the 
photon. From Eq. (|A5|) we have Mb = M a + Mb, and (we do not consider the polarization 
here) the divergent part of Z) sp i ns |M#| 2 is contained in the first two terms of the square 
brackets of 

2 



E \ M « 



2 e 2 G 2 v 2Mi x^fl-e pi-e 



spins 



2 Mjmm 



E 



V e 



I ■ k pi ■ k 



x 



D x EEl - D 2 \ ■ (p 2 + 1) - D x Euo - D 2 l ■ k 



(C2) 



Because 1 • (p 2 + 1) = lp2yo + 1 2 — r]/2 (choosing the z-axis along 1), we can split dT B as 



dT 



1 G 2 V dE dE 2 dVL t 
u ' 2~2 (2vr) 8 
21 - Pl 



a ( ,1 

4Mi- < A - lim / dn 

m 2 M 2 



1 d 3 kd 3 p v 
2n uj E v 



^ (Pi~P2-l-Pu- k) 



( Pl -k)(l-k) (I ■ k) w( Pl -k)\ 
1 A-x l d 3 kd 3 p v 4 fl-e Pl -e s 

4 Jo 2n uj E v £ \l ■ k pi ■ k 



x 



-DyE.c - D lk • -77 



^[|M b | 2 + 2Re[M a ] [M 6 ; 



(C3) 



where A x is given in Eq. (|B2|). In the first integral we perform the sum over polarizations 
indicated in Eq. ( p2| ) in covariant form. This allow us to identify this integral with the 
divergent integral I of Ref. U , namely, 



1 fflmax r 

I (E,E 2 ) = -Mm / dry 2(/ • Pl )I n (l, Pl ) - m 2 I 20 (l, Pl ) - M 2 I 02 (l, Pl ) 

4 J\2 L 



with 



(l,Pi) 



1 f d 3 p v d 3 k S 4 (pi — P 2 — I — Pv — k) 
2^ 



E u uj (l-k) m ( Pl -k) 
which are invariant integrals. The final form of Iq is 



J (E, E 2 



arctanh/? 



mrf m 



21n(^)+ln( 



4r_ 



, 777. 

2h,l T 



hi 



2mE® (q 2 — m 2 ) 



Ar 4 



(C4) 



(C5) 



(C6) 
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where 



(E + l)r±= [E°J 2 (q 



± 



Ell 2 (q 



m 

2 



a 2 E/A 



m 



a E/A - mV/16 



1/2 



(C7) 



and 



q 



Ml 



(C9) 



For the FBR of DP the lower limit of the dT B of Eq. flClp is r] min = 2lp 2 (yo — 1), and 
then Iq is no more infrared- divergent. We can calculate it from Eq. (|CJ4j) by changing A 2 by 
^mm- The second part of dTs in Eq. (C3) is infrared-convergent. It is more convenient to 
transform it to the form of Eq. (38) of Ref. f3[. We can accomplish this by performing the 
5-integration and by changing the rj variable to the y variable using 77 = 2p 2 l(yo — y)- We 



do not reproduce the result here because it is very long. It is given in our Eq. 



for dT 



b- 
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TABLE I. Values of the form factors used in our numerical calculations. For the first three 
decays we take the form factor ratios of Ref. [6] while for the other decays we use Ref. [13]. For 
A+ -» A we take the values of Ref. [10]. For convenience we include in the last column the 
uncorrected a° of Eq. (77) corresponding to this choice of form factors. 



Decay 


h 


h 


9i 


a 


n -> P 


1.000 


1.970 


1.261 


-0.0850 


A -> P 


1.236 


1.199 


0.890 


0.0200 


E~ — > n 


1.000 


-0.970 


-0.340 


-0.6319 


S- ->■ A 


0.000 


1.172 


0.601 


-0.7030 


£+ ->■ A 


0.000 


1.172 


0.601 


-0.6474 




1.225 


-0.074 


0.354 


0.2579 




0.707 


1.310 


0.899 


-0.1989 


3°^£+ 


1.000 


1.853 


1.267 


-0.1913 


A+^A 


0.350 


0.090 


0.610 


-0.9513 



TABLE II. Percentage radiative corrections [that is, Eq. (83) multiplied by 100] of the 
spin-asymmetry coefficient of the charged lepton in hyperon semileptonic decays. The predic- 
tion in the fourth column for A+ — > A uses the approach of Ref. [13], but it was not actually given 



there. 


Decay 


5a R = af — «q 


5a T = a J — «o 


5a Ref. [13] 


5a Ref. [6] 


n^P 


0.0119 


0.0095 


0.0101 




A^ P 


0.0813 


0.0014 


-0.0023 


-0.0 


£ - -► n 


0.0832 


0.0815 


0.0758 


0.1 


£- -»• A 


0.1432 


0.0836 


0.0770 




S+ -> A 


0.1287 


0.0755 


0.0911 






0.1024 


-0.0246 


-0.0310 




3- 


0.3327 


0.0371 


0.0212 




3° S+ 


0.3312 


0.0350 


0.0208 




-»• A 


0.0757 


0.1294 


0.1098 
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